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A new calculation of the neutron electric form factor is described. Preliminary 
results for the connected part on a 20 3 X 32 lattice are presented. The methods 
for calculating the disconnected part are briefly described. 

1 Introduction 

Previous lattice calculations of nucleon electromagnetic form factors on 16 3 x 24 
and 24 x 12 x 12 x 24 lattices 1,2 have had decent signals for the proton electric 
and magnetic form facors, G p e (q 2 ) and G P n {q 2 ) ) as well as the neutron magnetic 
form factor, G^(</ 2 ), but very bad signals for the neutron electric form factor, 
G"(q 2 ). The most that can be said is that the lattice results agree in sign and 
magnitude with the expected experimental results. If we are to move beyond 
the qualitative stage in comparison of lattice data to experiment, this situa- 
tion must be improved, especially considering the new CEBAF experiments 
planned 3 . A new, more precise, lattice calculation of G"(q 2 ) is described here. 
The lattice used is larger (20 3 x 32), more configurations (50) are employed 
and better analysis methods are used. 

2 Simulation Details 

Using the usual point nucleon interpolating fields 1 , one needs both the two and 
three-point Green's functions: 

G nn (t;p,T) = J2^ ¥ ^a-a{y^\T(xa(x)x a -(0)) |vac), (1) 

X 

Gnj^ih^qX) = -i^e-^ i r«,«(vac|T(x«(*2)J M (*i)x«'(0))|vac)(2) 

S 2 ,Si 

One now forms the connected and disconnected Wick contractions of the proton 
three-point function. The long Euclidean time limits of these finctions are 

fi (f-nVA E + m \Z\ 2 a 6 Et ( o\ 
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G n j in (t 2 ,t i; ir 4 ) ^ ^e-^-^e- E ^ (M±m) G n {q 2y (4) 

We do not measure these correlation functions directly, but instead analyze 
the ratio 

C n ( 2 )= ( 2£ ( G n j in (t 2 ,ti;q,T 4 ) \ ( Gnn^Uy y 

> {E + mJ \ G nn (t 2 ;0,T 4 ) J \G nn (t; q, T 4 ) J ' [ °> 

for G". Fixing t and t 2 , we look for regions in t\ where the signal stands out 
from the noise. (The 4x4 matrix r 4 = | ( □ o )•) 

Some results for the connected part of the neutron electric form factor are 
presented in Fig.l at k = .152 (« cr = .1564(2)). In the following we are at 
f3 = 6.0 working on 50 configurations, thermalized by 5000 pseudo-heatbath 
sweeps and separated by 1000 sweeps. The neutron source is at lattice time site 
6 and the zero- momentum sink is at time site 21. The solid and dashed lines in 
Fig.l are the "dipole" and "Galster" parameterizations of G"(q 2 ), respectively, 
evaluated with the nucleon mass and magnetic moment measured at k = .152. 
The squares represent the previous results for the lowest four nonzero lattice 
momentums from Ref.l; the diamonds are the new results on the larger lattice. 
One can see that the error bars are much reduced. Surprisingly, the results are 
much lower than expected. 

One can easily show that the expectation value of the exactly conserved lat- 
tice current density, J ^(x) = (J(x),ip(x)) = in[ip(x + a |[1 )(l + ■y II )Uj l (x)ip(x) — 
ip(x)(l — ■j ll )U ll {x)'il>(x + dp)], has a vanishing imaginary part, configuration 
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by configuration. Thus, in forming the disconnected amplitude one is looking 
at the statistical correlation between the imaginary part of the charge density 
and the imaginary part of the nucleon two point function, both of which are 
simply noise by themselves. 

Because of the average over gauge configurations, we may replace the point- 
to-point propagator M -1 (i"', t' ; x, t) with a version which is summed over all 
inital space-time positions in gauge invariant quantities as a consequence of 
Elitzur's theorem. This method 4 (the "volume method") allows the global 
measurement of the electromagnetic current, ^2g Jfi(%,t), necessary for the 
disconnected part. 

The other useful technique for disconnected amplitudes is the Z(N) noise 
method 5 . Numerically, I have found that Z(N) noise works better than the 
volume method (the same as SU (3) noise) for local operators such as ip(x)ip(x), 
but the volume method is better for the nonlocal J^x). Since interest here is 
directed at the disconnected electromagnetic contribution to the nucleon, the 
volume method has been chosen for the simulations. 

3 Conclusions, Remarks and Acknowledgements 

The method of calculating nucleon form factors on the lattice has been reviewed 
and a new calculation has been described. An excellent signal is produced 
for the connected G"(q 2 ) with error bars significantly smaller than previously. 
Although the chiral extrapolation has not been done, one can already conclude 
that the connected results are much smaller than expected. 

The method used for the disconnected calculation has been explained. One 
must look at the correlation of the imaginary parts of the nucleon propagator 
and charge density. The signal will be hard to detect, but must be present if 
the phenomenological results are to be recovered from QCD. 

This research is partially supported by NSF Grant No. PHY-9401068 and 
NCSA. I thank Phillip Kalmanson for helping verify the results presented here. 
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